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We consider quantum entanglement between gauge fields in some region of space A and its comple- 
ment B. It is argued that the Hilbert space of physical states of gauge theories cannot be decomposed 
into a direct product Ha (SiTCb of Hilbert spaces of states localized in A and B. The reason is that 
elementary excitations in gauge theories - electric strings - are associated with closed loops rather 
than points in space, and there are closed loops which belong both to A and B. Direct product 
structure and hence the reduction procedure with respect to the fields in B can only be defined if 
the Hilbert space of physical states is extended by including the states of electric strings which can 
open on the boundary of A. The positions of string endpoints on this boundary are the additional 
degrees of freedom which also contribute to the entanglement entropy. We explicitly demonstrate 
this for the three-dimensional Z2 lattice gauge theory both numerically and using a simple trial 
ground state wave function. The entanglement entropy appears to be saturated almost completely 
by the entropy of string endpoints, thus reminding of a "holographic principle" in quantum gravity 
and AdS/CFT correspondence. 

PACS numbers: 11.15.Ha; 03.65.Ud; 89.70.Cf 



Geometric entanglement entropy [1, 2, 3] of confining 
gauge theories has recently become a subject of extensive 
studies, mostly due to the discovery of its non-analytical 
behavior with respect to the size of the region. This 
non-analyticity was first predicted in the framework of 
AdS/CFT correspondence in [4, 5], and then within the 
Migdal-Kadanoff approximation in lattice gauge theories 
[(>]. Finally, a signature of nonanalytic behavior of en- 
tanglement entropy was found in numerical simulations 
of SU (2) lattice gauge theory [7]. 

In order to define the geometric entanglement entropy 
of some region A, one should represent the Hilbert space 
H of the field theory under consideration as a direct prod- 
uct of Hilbert spaces Ha and Hb of states localized in A 
and its complement B [8]: H = Ha ®Hb- For instance, 
for scalar field theory the spaces Ha and Hb can be con- 
structed by acting on some proper initial state with field 
operators localized within A and B. In other words, el- 
ementary excitations of scalar fields are associated with 
points in space, and any point can be classified as be- 
longing either to A or B (except for the set of points on 
the boundary, which has zero measure). One then de- 
fines the reduced density matrix for the fields in A as a 
partial trace of the density matrix of the ground state 
0) of the theory over Hb- Pa = Tr b |0)(0| . Geometric 
entanglement entropy S [A] is the von Neumann entropy 
of the reduced density matrix pa [8]: 

S[A] = -Tr A p A lnp A (1) 

However, in none of the works [4, 5, 6, 7] such a direct 
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product structure was explicitly constructed for gauge 
theories or even used. The papers [4, 5] mainly rely on the 
conjecture of [9] which relates the entanglement entropy 
of conformal field theories living on the boundary of AdS 
space with certain minimal hypersurfaces in the bulk of 
this space. The results of [G, 7] are based on the relation 
between entanglement entropy and the free energies of 
the theory on the spaces with multiple cuts, which was 
originally derived in [2, 3] for scalar field theories. 

The aim of the present work is to explicitly construct 
a decomposition of the Hilbert space of states of gauge 
theories into a direct product Ha®Hb- One of our main 
assertions is that such a decomposition is impossible for 
a Hilbert space of physical states, i.e. a Hilbert space of 
states which satisfy the Gauss law. The reason is that 
in contrast to scalar field theories, in pure gauge theories 
elementary excitations are not associated with points in 
space, but rather with closed loops, which arc the lines of 
electric flux or electric strings [10, 11, 12]. Correspond- 
ingly, closed loops cannot be classified as belonging either 
to A or B, since the set of loops which belong to both A 
and B has nonzero measure in the space of all loops. 

For the sake of brevity let us consider the simplest case 
of Z2 lattice gauge theory in (2 + 1) space-time dimen- 
sions. The arguments presented below can be generalized 
to other gauge theories in a straightforward way. The 
Hamiltonian of this theory is: 

£ = ^rI>-&En> ( 2 ) 

' 9 1 p iep 

where (3 g is the coupling constant, ^ and ^ denote sum- 

l P 

mation over all lattice links and plaquettes and zi and 
a 1 are the operators with eigenvalues equal to ±1 which 
satisfy {07, 1;} = 0. The Hilbert space is a space of all 
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FIG. 1: Basis quantum states for the simplest two- 
dimensional lattice. Thick lines are the links with mi = 1. 



functions [z{\ of ./^-valued link variables zi on a two- 
dimensional lattice. A convenient basis in this space can 
be parameterized by a set of labels mi — 0, 1 defined 

on each link of the lattice: ip[ z i',mi} ~ II ^z™'- This 

I 

basis is orthonormal with respect to the scalar product 
(01 102) = 01 [zi] tp2 [zi]- The Hilbert space of phys- 
{2i} . 

ical states Ho is obtained by imposing the Gauss law 
as a first-class constraint, which leaves only the wave 
functions which are invariant under gauge transforma- 
tions: [zp] = [zi]. The basis states in Ho can be la- 
belled by such sets {mi}, for which the links with mi — 1 
form closed non-self-intersecting loops, and can be ob- 
tained by acting with all possible Wilson loop opera- 
tors W [C] = Yl z i on t ne trivial strong-coupling ground 

state with [z{\ = const. The operator W [C] thus cre- 
ates an electric string on the loop C. 

As the simplest example consider the lattice which con- 
sists of just two plaquettes, p\ and P2 (see Fig. 1). The 
Hilbert space of physical states on such lattice is four- 
dimensional and contains the states with electric flux 
flowing cither around plaquettes p\ or P2, or around both: 

-01 [Zl] = 1, 02 [Zl] = Z pi , 03 [zi\ = Zp 2 , 04 [zi] = Z pi Z p21 

where z p = Yl z i- By direct counting of states one 

can make sure that the Hilbert space spanned on this 
four states can not be represented as a direct product of 
Hilbert spaces of functions on any complementary sets of 
links. 

A proof of this assertion in the general case is also 
rather obvious. Let us denote the Hilbert spaces of all 
functions on links within A or B as Ha and Hb- We 
would like to decompose the Hilbert space of physical 
states Ho as Ho — Ha®T~Cb, so that Ha and Hb are the 
subspaces of Ha and Hb- Consider some state with the 
wave function of the form [zi] ~ Y\ Zi, where the C is 

c 

some closed loop which belongs both to A and B. This 
state is a direct product of two states in Ha and Hb with 
wave functions ipA [z{\ ~ n z i and ips [zi] ~ n z h where 

Ca Cb 

Ca and Cb are the parts of the loop C which belong to A 
and B. Since the basis states corresponding to different 
loops are orthogonal, this decomposition is unique, and 



the states ipA and 0s should necessarily belong to Ha 
and Hb- The spaces Ho and hence Ha and Hb should 
also contain the trivial strong-coupling ground state with 
0o [zi] = const = 0AO0BO, 0a,bo [zi] = const. Therefore 
if Ho is indeed a direct product of Ha and Hb , it should 
also contain a direct product of, say, 0a with ipB- The 
wave function for such a state is proportional to Y[ z i-> 

Cb 

but Cb is in general not closed, therefore the Gauss law 
is violated at the boundary between A and B. We thus 
arrive at the contradiction, which completes the proof. 

Let us now try to extend the Hilbert space of physi- 
cal states in some minimal way, so that a direct product 
structure can be introduced. Consider again the lattice 
on Fig. 1. Since the cause of the trouble are the strings 
which cross the boundary between A and B (this is the 
link in the middle on Fig. 1), let us try to cut such strings 
along this boundary. This yields four more states, which 
are schematically shown on Fig. 1, b. The resulting ex- 
tended Hilbert space is now eight-dimensional. Assume 
that the link in the middle and all the links to the right 
of it belong to the region A, and that all other links be- 
long to the region B. It is easy to see that the extended 
Hilbert space can be represented now as a direct prod- 
uct of a two-dimensional Hilbert space of functions on 
the links in B times a four-dimensional Hilbert space of 
functions on the links in A, which are not invariant under 
the gauge transformations only at the endpoints of the 
link in the middle. In the general case, the minimal ex- 
tension Ho of Ho is also the Hilbert space of states which 
violate the Gauss law only at the boundary OA between 
A and B. Minimal extension here means that the Gauss 
law does not hold in a minimal number of points, and 
Ho is a minimal subspace of H which contains Ho and 
which can be represented as a direct product. To see that 
this is indeed so, suppose that the Gauss law still holds in 
some points on dA. In this case one can repeat all the ar- 
guments above for the loop C which goes through these 
points and arrive at the same contradiction. It should 
be mentioned that such an extension of Hilbert space of 
physical states was already implicitly used in the study of 
topological entanglement entropy [13, 14] and entangle- 
ment entropy of spin networks [15], although only for a 
particular class of ground states which can be thought of 
as the ground states of lattice gauge theories in the limit 
of infinitely strong coupling. In the extended Hilbert 
space one can naturally define the partial trace over the 
fields in B as a sum over all links which belong to B: 

(z' l \Tr B p\z' l ') = p[M},MY\ (3) 

{zi},ieB 

If the expression (3) was simply postulated without intro- 
ducing the extended Hilbert space, it would lead to pe- 
culiar density matrices which satisfy the Gauss law, i.e., 
which commute with the generators of gauge transfor- 
mations, but which nevertheless can not be represented 
in the form pa = J2 w i IV'iXV'il j where \ipi) are physical 

i 

states. As an example, consider again the state with the 
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FIG. 2: Thin tubes of electric flux in the vicinity of a charged 
black hole with charge Q = 4. 

wave function ip [zi] ~ Y[ z i > where C is some loop which 
lec 

belongs both to A and B. The reduced density matrix 
Pa [z'i, z['] ~ [ is invariant under simultaneous 

l£C A 

gauge transformations of z[ and z'{ , but obviously can 
not be decomposed over the basis of physical states. On 
the other hand, in the extended Hilbcrt space such de- 
composition is trivial. Thus one could in principle arrive 
at the concept of the extended Hilbcrt space by trying to 
interpret the outcome of the reduction procedure (3). 

We thus conclude that in order to define the spaces of 
states of gauge fields which arc localized cither in A or B 
one should extend the Hilbcrt space of physical states and 
include also the states of electric strings which can open 
on the boundary dA between A and B. This boundary 
can be therefore considered as a sort of D-brane for elec- 
tric strings, with the positions of their endpoints being 
the additional degrees of freedom which emerge as a re- 
sult of such extension. It is reasonable to conjecture that 
these extra degrees of freedom should contribute some- 
how to the entropy of entanglement between A and B. 
In order to understand how this can happen, it is in- 
structive to consider a confining Abelian gauge theory, 
such as compact QED in the strong coupling limit [10], 
which interacts with a charged black hole with charge 
Q (see Fig. 2). Since the region inside the horizon is 
not accessible to an external observer, in order to find 
the density matrix of the fields outside of the black hole 
one should trace over all states of fields inside the hori- 
zon, which is one possible way to relate the Bekenstein- 
Hawking entropy with the entanglement entropy [16, 17]. 
For confining theory the charge should be connected with 
the horizon with \Q\ thin tubes of electric field, and \Q\ 
such tubes should also go out of a black hole. The po- 
sitions of the endpoints of these flux tubes outside and 
inside of the horizon should be completely uncorrelatcd, 
otherwise such a system would violate the "no - hair" 
theorem. On the other hand, before the formation of 
a black hole all flux tubes were continuous, and ended 
only on the charge Q. The entropy associated with the 
uncertainty in the positions of endpoints on the horizon 
should therefore also contribute to the entropy of a black 
hole and hence to the entanglement entropy of the fields 



behind the horizon. 

In order to estimate the contribution of these addi- 
tional degrees of freedom on the boundary dA, let us 
first calculate the entanglement entropy for a simple trial 
ground state wave function of Z 2 lattice gauge theory in 
(2 + 1) space-time dimensions. This trial ground state 
is a superposition of all possible configurations of closed 

electric strings with the weight exp ^— ^ ^2m^j, where 
m i is the total length of strings: 

1 

*oN = C ex P( ~f ^2mi\^[zi;mi] (4) 

where Smi = means that we sum only over closed elec- 
tric strings and C is the normalization constant. For- 
tunately, the trial ground state (4) yields an expression 
for the entanglement entropy which has a universal form, 
and which, as wc will see from the results of numerical 
simulations, seems to be a good approximation for the 
entanglement entropy of a true ground state of Z 2 lat- 
tice gauge theory for all values of the coupling constant. 
The wave function (4) would be a product of functions 
of individual link variables if the constraint Smi = was 
omitted. Thus in some sense the trial ground state (4) is 
only "minimally" entangled because the electric strings 
are closed, and there is no entanglement between closed 
string states in the bulk of the regions A and B. Despite 
its simplicity, the wave function (4) can qualitatively de- 
scribe the most essential features of the dynamics of elec- 
tric strings in Z 2 lattice gauge theory, such as the emer- 
gence of percolating strings in the vicinity of a quantum 
phase transition, due to the dimensional reduction prop- 
erty [18]. 

The entanglement entropy for the state (4) can be cal- 
culated using the "replica trick" S [A] = — lim Tr p B A 

s— »1 

[2, 3]. Applying the reduction procedure (3) and using 
the orthogonality of the basis states ^[zi;m{\, one can 
rewrite the trace Tr p s A as: 

Trp A = C 2s J2 ■■ E cxp( -a ( 5 ) 

Here one should sum only over such {mik} for which 
the electric strings in A described by {mik}, I S A and 
{mik+i}, I € A end in the same points on dA as the 
string in B described by {mik}, I E B. Iterating this 
constraint for all k = 1, . . . , s, we see that we actually 
sum over s independent string configurations in A and 
B, which all intersect dA in the same points x±, . . . , x m . 
It is convenient to sum separately over the coordinates 
of these endpoints and over all configurations of electric 
strings in A and B which end in x±, . . . , x m . After such 
rearrangement the sums over {mik}, I S A and {mik}, 
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I e B in (5) factorize: 

{xi,...,x m } \{mi} \ ( 



(6) 



where the prime over the sum in the brackets means that 
that the electric strings described by {mi} should inter- 
sect OA only in the points {x\, . . . , x m }. Remember now 



that exp y—a E m i J is the square of the trial wave func- 
tion (4) in the basis of states ip [zi; mi], or, in other words, 
simply the probability distribution of the configurations 
of electric strings on the lattice. The constrained sum 
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i 



is then nothing but the prob- 



{nn} 

ability distribution p [{x\, . . . , x m }] of the intersection 
points between electric strings and the boundary dA in 
the ground state of the theory. Taking now the derivative 
of (6) at s = 1 and using the replica trick, we find that 
for the trial ground state (4) the entanglement entropy is 
simply the classical Shannon entropy of the probability 
distribution p [{x\ , . . . , x m }] : 



S[A} = - Yl pK*i.- 



J] lnp [{si,. 



• }](7) 



i } 



Although this result was derived only for a particular 
wave function (4), its simple and universal form suggests 
that it can be applied to any quantum state, for which 
the entanglement between the states of the closed strings 
in the bulk of the regions A and B can be neglected. An 
expression similar to (7) has been also obtained recently 
for the entanglement entropy in loop quantum gravity 
[15]. 

For the purposes of numerical simulations, it is con- 
venient to represent the expression (7) in a somewhat 
different form, assuming for simplicity that for a given 
m all the intersection points xi, . . . ,x m are not corre- 
lated. Taking into account that only one electric string 
can pass through one lattice site, after some simple trans- 
formations one can rewrite the expression (7) as: 



S [A] = So [P(m)} + (In 



\dA\\ 



ml {\dA\-m)\ 



) (8) 



where So [P (m)] is the entropy of the probability distri- 
bution P (m) of the number of intersection points and 
\dA\ is the total length of the boundary dA in lattice 
units. 

In order to check the validity of the expression (7) 
for the true ground state of Z2 lattice gauge theory, we 
have measured the entanglement entropy and the en- 
tropy of the probability distribution of intersection points 
p \{x\, . . . , Xm}] in lattice Monte-Carlo simulations. The 
region A was a square of size ! X L To calculate the 
entanglement entropy, we have used the same method 
and the same approximations as in the work [7], mea- 
suring only the differences of entropies S [A 1 ] — S [A] = 




FIG. 3: Gluing of s = 3 (1 + 1) -dimensional square lattices 
into a three-sheeted Riemmann surface. 
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FIG. 4: Lattice derivatives of the entanglement entropy and 
the entropy of intersection points over the size of the region 
A at p g = 0.788, p s = 0.21, 16 3 lattice. 



5(7 + 1) — S (I). Let us note here that the arguments 
presented in this Letter actually justify the method used 
in [6, 7] to calculate the entanglement entropy of lat- 
tice gauge theories. The lattices with cuts used in [6, 7 
have the topology C( s ) ® J D ~ 2 [2, 3], where C( s ) is 
the s-sheeted Ricmann surface, T D ~ 2 is the (D — 2)- 
dimcnsional torus and D is the dimensionality of space- 
time, and were constructed by gluing together s square 
lattices along the links which belong to the region A (see 
Fig. 3) . Such construction essentially relies on the notion 
of the partial trace (3) over the links in B , which only 
makes sense for the extended Hilbert space constructed 
above (see [2, 3] and Appendix A in [7] for a detailed 
derivation) . 

Similarly to the entanglement entropy, the positions 
of intersection points between electric strings and the 
boundary dA are not directly measurable in Monte-Carlo 
lattice simulations. The reason is that the correspond- 
ing operator is not diagonal in the eigenbasis of the link 
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FIG. 5: Lattice derivatives of the entanglement entropy and 
the entropy of intersection points over the size of the region 
at I — 8 at different values of the coupling constant f3 g for 16 3 
lattice. 



variable operator zi. Fortunately, for Zi lattice gauge 
theory in (2+1) dimensions one can use the Kramers- 
Wannier duality transformation and express the expec- 
tation values of field operators which are diagonal in the 
basis ip [zi]m{\ as expectation values in the dual theory 
[10, 12], which is the Ising model in (2 + 1) dimensions 

with the Hamiltonian H = ^i- ^ ^s~Ps X) z s z s+ i, where 

s i 

^2 denotes the sum over all lattice sites. The configura- 

s 

tions of electric strings in the gauge theory can be identi- 
fied with the closed loops formed by the links which are 
dual to the links with z s z s+ i = —1 in the Ising model, if 
the coupling constant (3 S is related to the coupling con- 
stant (i g of the gauge theory as th/3 s = exp(— 2(3 g ) [12]. 
Thus in order to find the entropy of the probability dis- 
tribution of intersection points between electric strings 
and dA we have performed Monte-Carlo simulations of 
the three-dimensional Ising model at the dual coupling, 
extracting the configurations of electric strings from the 
configurations of spin variables z s on some fixed time 
slice as described above. Since it is extremely difficult 
to measure exactly the m-point probability distributions 
p m (xi, . . . ,x m ), we have neglected the correlations be- 
tween xi, . . . , x m and used the expression (8). 

Lattice derivatives S (I + 1) — S (I) of the entanglement 
entropy and the entropy of the distribution of intersec- 
tion points are plotted on Fig. 4 and Fig. 5. For the plot 
on Fig. 4 the coupling constants /3 g = 0.788, (3 S = 0.21 
are close to the critical point with (3 SC « 0.22, at which 
percolating electric strings emerge. It is reasonable to ex- 
pect that for percolating strings the intersection points 
{x\, . . . , x rn } are indeed uncorrelated. At these values 
of coupling constants the differences S (I + 1) — S (I) of 
both entropies indeed agree within the error range and 
practically do not depend on /, which means that the en- 
tropy is proportional to the area of the boundary with a 
good precision. We thus recover the familiar "area law" 



for the entanglement entropy [19]. The data plotted on 
Fig. 5 was obtained for the fixed size of the region A 
I = 8 and for different values of the coupling constants 
(3 g and (3 S . Again, only the differences S (I + 1) — S (I) 
were compared. The dependence on j3 g is similar for both 
entropies, however, in the strong coupling phase of the 
lattice gauge theory the entropy of intersection points is 
somewhat larger than the entanglement entropy. This 
discrepancy can be explained cither by correlations be- 
tween different intersection points or by some systematic 
error of the approximation used to calculate the entan- 
glement entropy [7]. It is not likely that the difference 
between the two entropies can be caused by additional 
entanglement between lattice degrees of freedom in the 
bulk of the regions A and B, since in this case it is rea- 
sonable to expect that the entanglement entropy should 
be larger than the entropy on the boundary. 

It should be stressed that our results does not exclude 
the entanglement of states of electric strings which do 
not cross the boundary between A and B, however, in 
this case the pattern of entanglement is the same as for 
the scalar field theories in the bulk of A and B. Such 
situation may be realized when the entanglement en- 
tropy of gauge theories is calculated in the framework 
of AdS/CFT correspondence using the prescription of 
[4, 5, 9]. In the quasiclassical approximation, in which 
the entanglement entropy is proportional to the minimal 
area of some hypersurface, the full string theory on AdS 
space is reduced to supergravity, which is the theory of 
particles rather than strings. Since the strings in AdS are 
believed to describe electric strings in the dual gauge the- 
ory, this means that the latter are very small and we are 
in the strong coupling limit (cf. the strong coupling limit 
of lattice gauge theory, where the dominant excitations 
are the electric strings which wind around one plaquette 
only). Correspondingly, there are very few strings that 
cross the boundary between A and B, and the "brane" 
effect described above should be negligibly small. 

We thus see that the expression (7) works well for all 
sizes of the region A and for all values of the coupling 
constant, and the entanglement between the states in 
the bulk of the regions A and B seems to be rather small. 
The situation indeed looks as if the entanglement entropy 
was saturated by the uncertainty in the positions of the 
endpoints of open electric strings in the extended Hilbcrt 
space at the outer side of dA. In view of the actively dis- 
cussed "holographic descriptions" of field theories, which 
relate the properties of string theories in the bulk of AdS 
space or modifications thereof with the properties of the 
field theories on its boundary [20, 21, 22], it seems rather 
tempting to associate the entanglement entropy of the 
region A with the classical entropy of some statistical 
theory on its boundary. The endpoints of electric strings 
could then be interpreted as elementary degrees of free- 
dom in such a theory. This can have interesting conse- 
quences in black hole physics - for instance, a confining 
gauge theory in the vicinity of a black hole might induce 
some two-dimensional field theory on its horizon, simi- 
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larly to the scenario discussed in [23]. It could be also 
interesting to generalize the expression (7) to the case of 
non-Abelian gauge theories, which is not a straightfor- 
ward task, since the notion of conserved electric flux is 
not well understood for non-Abelian gauge groups. 
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